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I. INTRODUCTION 



The development of quantum theory of Brownian motion laid the foundation of non- 
equilibrium statistical mechanics of a system coupled to its environment [U [2]. This has 
significantly affected the course of research in quantum optics, condensed matter physics 
and chemical dynamics, over the last several decades[3j 0]. Although its major thrust lies 
on calculation of non-equilibrium properties in terms of relevant correlation functions, a 
few interesting attempts have been made in evaluating thermodynamic quantities of small 
quantum system coupled to its environment which exerts a finite dissipation on the system. 
For example, the free energy of a quantum oscillator interacting via dipole coupling to black- 
body radiation field modes has been calculated [5] to obtain a temperature-dependent shift 
in the free energy resulting in modification of Planck's formula. The effects of dissipation on 
thermodynamic functions have been explored for quantum oscillators in contact with Ohmic 
and non-Ohmic reservoirs[6] and for the magnetic moment of an electron gas[7]. It has also 
been shown that while a classical free particle does not obey third law of thermodynamics jH] 
its coupling to a thermal reservoir renders quantum nature and allows recovery of third 
law[5H8]. That quantum dissipation helps to ensure the validity of third law has been a 
major finding for a dissipative quantum oscillator jU], for which the specific heat at low 
temperature exhibits power law behaviour. 

The focus of the present paper is the calculation of thermodynamic functions, partic- 
ularly the specific heat of a spin-| particle coupled to a spin bath. To bring forth the 
discussion in an appropriate perspective we first note the following points. When a small 
system is made open by coupling it to an environment the system experiences a dissipative 
flow of energy into the environment. It is easy to anticipate that such a flux would give 
rise to a dynamical correction over the usual thermodynamic contribution responsible for 
the temperature dependence of the system, a situation reminiscent of Kramer's theory of 
activated rate processes[TJ [2]. In the later case the rate constant is essentially a product of 
the equilibrium contribution due to transition state theory and a dynamical factor due to 
the frictional coefficient characteristic of the thermal bath [I]. It is possible to realize the 
nature of dynamical corrections over the thermodynamic functions by spatially restricting 
the system size when the fluctuations of the internal energy of the system can be quite 
significant. In this paper we focus on the dynamics of an open spin system experiencing 

2 



spin-bath induced dissipative effects [9] as an underlying microscopic description of this be- 
haviour. To be specific, the behaviour of specific heat for small systems can be analysed 
by examining the finite size effects on thin films and fine particles [T0HT3] . The strategy for 
calculation is based on the free energy of the spin system coupled to the spin bath [TU EE] 
minus the free energy of the spin bath in absence of the system. The desired free energy 
of the interacting spin system turns out to be an integral over the free energy of a single 
system multiplied by a density of states related to the susceptibility (16] derived explicitly 
from the associated c-number quantum Langevin equation. For implementing the scheme 
we have employed the spin coherent state representation [T7] of our proposed Hamiltonian 
of the spin-spin-bath model to cast the c-number Hamiltonian into an oscillator-oscillator 
bath model. The key point is the use of Holstein-Primakoff transformation [18] which sets 
up a mapping between spin and boson operators. Interestingly this allows us to recover the 
traditional form of quantum Langevin equation (its noise correlation, however, is guided 
by the spin characteristics of the bath) and consequently the associated susceptibility. We 
derive the explicit analytic expression of the specific heat and its high and low temperature 
limits, which carry the signature of dissipative effect of the environment. It has been shown 
that specific heat can be expressed as a sum of the products of an equilibration factor that 
contains the usual temperature dependence and a dynamical correction factor characteristic 
of the dissipative flow of energy from the system. Finally, our analysis is correlated with 
experiments dealing with finite size effects on magnetic ordering temperature in the specific 
heat measurements of CoO/Si0 2 monolayers [TU]. 



The outline of the paper is as follows: We introduce the spin-spin-bath model and the 
associated Hamiltonian followed by its spin coherent state representation in Sec II. The basis 
of this analysis is the Holestein-Primakopff transformation which further allows us to map 
the problem to an oscillator-oscillator-bath Hamiltonian. The c-number quantum Langevin 
equation is a direct consequence of this mapping. In Sec III we calculate the thermody- 
namic functions, free energy and specific heat for the interacting spin system and the high 
and low temperature limits of the derived expressions to correlate them with experimental 
measurement of specific heat on thin films. The paper is concluded in Sec IV. 
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II. QUANTUM DYNAMICS IN A SPIN BATH 



A. The Model 

To set up the problem of quantum dissipation of a two-level system in a sea of a two-level 
atoms, we consider the following Hamiltonian: 

H = hu S z + hJ2^kS zk + Hj29k\^i-(sl + S k }j (& + s)\ (2.1) 
k k ' 

Here the first term is the Hamiltonian for the system, which is specified by the Pauli pop- 
ulation difference operator S z , the other operators being S\ S and /, where and S are 
the usual raising and lowering operators and J is the identity operator. The second term 
corresponds to the reservoir Hamiltonian; the bath operators are denoted by S k , S k and S zk 
(the subscript refers to the k-th atom of the bath). The last term represents the interaction 
between the system and the bath atoms. u and u k are the characteristic frequencies of 
the two-level system and the k-th two-level bath atom. g k and g' k are the coupling constant 
and a frequency scale factor (this is a necessity for correct normalization of the spin coher- 
ent states used in the next section). The Pauli operators for the system follow the usual 
commutation relations as given by: 

[S\ S] = 2S Z ; [S\ S z ] = -S*; [S, S z ] = S (2.2) 

We have = S x + iS y , and S — S x — iS y . The commutation and anti-commutation rules 
between the spin-| operators are given by, 

SS j = 2S Z 

&S + S& = 1 (2.3) 

So, S z — N — |, where N is the number operator S^S. The anti-commutation relation in 
Eq. 2.3 has the immediate consequence that spin-^ particle or two-level atom obey Fermi- 
Dirac statistics. In view of Eqs.l and3 the system Hamiltonian H s may also be expressed 
as H s = Hu S z = hjjj (N — |), while the eigen- value equation for the number operator may 
be written as N\n) = n\n) with n=0,l. The general state with no quanta is denoted by 
|0) which satisfies JV|0) = and iJ s |0) = Hu (N - |)|0) = -^|0) and the state with one 
quantum is denoted by |1) which obeys N\l) = 1 and H s \l) = hu (N - = 
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Similar relations hold good for bath operators. Note that unlike harmonic oscillator, the 
Hamiltonian is asymmetric, both for the system and the bath and the Hilbert space is 
two- dimensional . 

Our object in the next section is to construct a c-number Hamiltonian from Eq.l using 
spin coherent state representation as introduced by Radcliffe[I7] a couple of decades ago. In 
what follows, we now give a brief outline of the coherent states of a general spin-S system, 
from which the spin-| coherent states can be obtained as a special limit. 

B. Spin coherent states 

We consider a single atom with spin S, characterized by the components of mutually non- 
commuting spin vector operators S x , S y and S z i.e. [S x , S y ] = iS z and cyclic permutations. 
The state of a spin is denoted by \S,m s ) which is the eigenstate of operators S 2 and S z 



so that we have S\S,m s ) = S(S + l)\S,m s ) and S z \S,m s ) = m s \S,m s ). The projection 
quantum number m s runs from -S, -S+1,...0,...S-1, S. We further define h = S — S z , such 
that the state \n) becomes an eigenstate. S\ S, n and S z have the following properties; 

S j \n) = {2Sfl^l- T ^^\n - 1) (2.4) 

S\n) = (2Sy/ 2 VT+^^l-^\n + 1) (2.5) 
h\n) = n\n) (2-6) 
S z \n) = (S-n)\n) (2.7) 



where n(= S — m s ) is the eigenvalue of n representing the difference between z-component 
of the spin and its maximum value, "n" is therefore allowed to run from to 2S. Now the 
state with minimal projection (for which m s = —S) is considered as a vacuum for a set 
of boson operators. Each subsequent state with higher projection quantum number is a 
boson excitation of the previous one. This can be achieved by considering the creation {a))/ 
annihilation (a) operator as, 

a)\ n ) = Vn + l\n + l) (2.8) 
a\ n ) = y/n\n - 1) (2.9) 
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By comparing Eqs (4-7) and Eqs (8-9) one obtains 



& = {2Sy'\ll-^a (2.10) 



S=(2S)^a^l-^ (2.11) 
where h = a) a (2-12) 

Eqs (10-12) are precisely the Holstein-Primakoff transformation [T8] which sets up a cor- 
respondence between the spin (-deviation) state and the bosonic state as: 



| - S + m s ) Spin i — > (a t )"|0) BosO n {m s = n) (2.13) 

Then the operator S = S x — iS y as given in Eq.ll creates spin deviations, when acts upon 
the newly defined ground state |0) (with minimal projection m s = —S). In fact we have, 

/ n !2<?l \ 1/2 

(S) n \0) = . ., \n) (0 < n < 2S) (2.14) 



(2S -n)\ 

Now with the help of Eq.14 the normalized coherent state is defined as: 

|A*> = N~ x ' 2 exp (jjJ?) |0) 

= (1 + \fi\ 2 )- s exp (ji£l\ |0) (2.15) 

The states defined by Eql5 do form a complete set[T7j. 

For the normalized coherent state given by Eql5, the following matrix elements may be 
derived [HI [19]: 

= (i + H 2 )" 2S H 2 ^(i + H 2 r 
-™K.= c(SM)\»\ 2 (2-16) 



i + \fi 

where C(S, = 

Similarly the other relevant matrix elements are given by 



1 + |/i 

and 



M#lA*> = T^TIa = C^HV ( 2 - 17 ) 



(^> = -^e C(SM)f (2-18) 
1 + \n\ z 
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We use the set of expressions (16-18) for calculation of matrix elements of the Hamiltonian 
(Eq.l) in the next section. Application of spin coherent states are well known in the context 
of ferromagnetic spin wave, phase transition in Dicke model of super-radiance, equilibrium 
statistical mechanics of radiation-matter interaction and so on. For details we refer to 
Klauder and SkagerstampO] . 



C. Spin-spin-bath Hamiltonian in coherent state representation 

Now we return to Eq.l and carry out quantum mechanical averaging with product sepa- 
rable coherent states of the system and the bath at t=0, \£,)\pi)\p2)---\hn), where |£) denotes 
the coherent state of the two- level system and l/i^) corresponds to the coherent state of the 
k-th bath atom. In using the coherent state representation of the Hamiltonian, we express 
the identity operator in Eq.l in terms of the system operators; 

I = S 2 + £ t2 + + &S (2.19) 

Here we have used the anti-commutation relation of the spin-^ or two-level system opera- 
tors. The consequence of expressing the identity operator asl9 will be clear immediately. 
The set of c-numbers for the bath degrees of freedom {/Zfc(O), a4(0)}, k = 1,2,.. .00, obtained 
from the different matrix elements, as discussed in Eqs. (16-18), can then be expressed in 
terms of momenta {/3k} and coordinates of bath oscillators, using the transforma- 

tion as /3k = z^ ^^ — (/4 — /ifc) and a>k = \J ~^^^ (a4 + /■**;)■ Similarly for the system 
oscillator co-ordinate r and momentum p can be expressed as r = v/ ^j^ (£* +0 and 
p = i\J (£* — £). Then the Hamiltonian in the coherent state variables can be ex- 

pressed in terms of these redefined quantities as, 

^4^ 2 +e!4i><--t) 2 (2 - 20) 

k k v K ' 

where Ck is given by 2^ gkg' k u)kU)od(\£\) or 2 J g (| |f| | ) g'k v^fc^o after discarding the irrelevant 
constant terms. The above Hamiltonian is different from our starting Hamiltonian oper- 
ator equationl because of the c-number nature of the coherent state variables. It is thus 
possible to reduce the quantum dynamics of a spin-| particle in a sea of spin-bath to the 
dynamics of a particle in an oscillator bath. This reduction is realizable in view of the well 
known kinship between the spin-| algebra and the algebra of bosons according to Schwinger's 



prescription [21]. Here, we use the Radcliffe coherent states, where Holstein-Primakoff trans- 
formation plays a connection between the spin-^ algebra and boson operators as does the 
Schwinger mapping|21J. Furthermore, high spin limit of a spin coherent state merges to 
the harmonic oscillator coherent state[17j. The spin-bath as a set of oscillators had been 
realized explicitly earlier [22] in a different context. The Hamiltonian20 is therefore the c- 
number equivalent of the Hamiltonian operator 1. Thermalization of the spin in a spin bath 
can thus be conveniently understood in terms of this c-number Hamiltonian. 

The equations of motion for the particle and for the bath variables according to the 
Hamiltonian Eq.20 take the form: 

Ck 



r + wj}r = ^2 °k y a k - ^2 r J ( 2 - 21 ) 

a k + bj\a k = c k r (2.22) 



Solving Eq.22 and substituting in Eq.21 with the condition r(0) = 0, we eliminate the bath 
degrees of freedom to obtain our desired c-number Langevin equation for the particle, 

r + [ dt'"/(t-t')r(t')+u%r = r](t) (2.23) 
Jo 

where, 

c 2 

7(t) = V -4 cosa; fc t (2.24) 



rj(t) = c k a k (0) cos uj k t + —/3 k (0)smuj k t (2.25) 
, u k 

k 

are the memory kernel and noise, respectively. 

In order to quantify the properties of the thermal bath, it is convenient to introduce a 
spectral density function J(uj) associated with the system-bath interaction. 

7T cl 



J( W ) = -V^( W - Wt ) (2.26) 

In terms of spectral density function J(u), one may rewrite the expressions for memory 
kernel in Eq.24 as, 

7 (t) = - f du^- cosut (2.27) 



while the noise rj(t) must satisfy the characteristics of the spin bath at equilibrium: 

(V(t)) s = (2.28) 

cosu(t-t') (2.29) 



■taring 



2 y 2KT' 



To ensure that the c-number noise rj(t) is zero-centred (Eq.28) and satisfies fluctuation- 
dissipation relation (Eq.29), it is necessary to assume a canonical distribution of Gaussian 
form for statistical averaging (...) s over c-number bath variables as follows: 



(2.30) 

This is the spin bath counterpart [TJJ US] of Wigner canonical thermal distribution function |23j 
for the harmonic or bosonic bath. Here N is the normalization constant. The width of dis- 
tribution is given by tanh(^^), which is related to the average thermal excitation number 
np(uJk) of the bath as tanh{^^) = 1 — 2np(u!k), np(u!k) being Fermi-Dirac distribution 



function. The width of Wigner distribution on the other hand is determined by coth( 2l r 
Note that at T=0, both distributions merge at a single value. The differences begin to 
appear at finite temperatures. At high temperatures, coth-factor reduces to a factor that 
results in the recovery of the classical limit. On the other hand the distribution (Eq.30) for 
the spin bath does not. Therefore the thermalization of the particle in a spin-bath can be 
described by this canonical thermal distribution30 and the c-number Hamiltonian (Eq.20). 



III. THERMAL EQUILIBRIUM OF A SPIN COUPLED TO A SPIN-BATH 

A. General expression for specific heat 

Here we derive the expression for specific heat (C v ) of a spin-spin-bath system using 
normal mode frequencies and examine the high and low temperature limits. We expect 
from the discussions in the previous section, that at low temperatures the spin-bath closely 
follows a bosonic bath while it differs at high temperatures. 

We begin by considering the normal mode solutions of Eqs.21 and22. To this end we 
write, 

r (t) = r {u)e- iult 
a k {t) = a k {u)e~ iwt k = 1,2, ... (3.1) 
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Using the solutions31 in Eqs21 and22 we obtain, 



(uj 2 - uj 2 ) r (u) = ^ c k ( 
k 



Ck , s 



UJ 



(ul - uj 2 ) a k (u) = c k r (uj) 
Elimination of ro(o>) from Eqs.32 and33 yields, 



, ,2 
k W k 



Id 



00* 



OJt 



(3.2) 
(3.3) 

(3.4) 



From the above equation the generalised susceptibility can be defined [24"] as, 



k(cj) 



{uj 2 - ujI) 



UJ 

9 S 
UJ Z — UJt 



-1 



(3.5) 



For discrete modes, k(u) has poles on the real axis at the normal mode frequencies of the 
system-plus-bath and has zeros at normal mode frequencies of the bath only. This allows 
one to write k(uj) in terms of the products of the ratios of the functions of the appropriate 



normal mode frequencies Tli(uj 2 — uj 2 )/Ylj(uj 2 



UJ 



i and j indices correspond to bath and 



system-plus-bath normal modes, respectively. Using Eqs.24 and26, the above equation can 
be rewritten; 

k(uj) = \—uj 2 + uj 2 — iuj^iuj)^ 1 (3.6) 

where j(cj) is Laplace transform of j(t). In what follows we show that k(uj) is responsible 
for dynamical corrections. We now consider that the spin coupled to the spin bath is in 
thermal equilibrium at temperature T. This system has well defined free energy F(T) which 
can be expressed [T6] as the difference between the free energy of spin-spin bath system and 
the free energy of the bath in absence of the spin. This can be related to the dynamical 
susceptibility k(uj) as follows; 



duf(uj, T)Im 



d 

duj 



where 



-KTln 



1 + exp 



Iii\k(u)\ 



Huj 
KT 



(3.7) 



(3.8) 



is the free energy of a spin-| particle or two-level atom at frequency uj. The thermodynamic 
functions, entropy and specific heat can be derived from the following relations: 

dF 



S 



(3.9) 



C V = T 



as; 
WT 
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— T 



d 2 F 
gf2 



(3.10) 



Now considering to heat bath to be Ohmic, i.e., 7(0;) = 70 one can simplify Eq.36 as 
follows, 

k(oj) — [— u 2 + Uq — icu7o] 1 (3-11) 



which results in 



Im 



(ln\n(oj)\) 



Then Eq.37 becomes, 



F(T) 



IT 



7Y 



duf(w,T) 



dujf(uj,T) 



70 (u 2 + up 
(w 2 - u 2 f + 7^2 



70 (u 2 + u 2 ) 

(u, 2 - co 2 ) 2 + 7o ^ 2 

z z* 

+ 



OJ 2 + Z 2 ' U 2 + 2* 2 



(3.12) 



(3.13) 



where, z = ^ + icui and 2* = ^ — iui and Wi 



2 



"8 - * ■ 



Eq.43 can be rewritten as, 



F(T) = - / dw/^T) ( + . . , 

n J \ u 2 + z 2 oj 2 + z* 2 



(3.14) 



Here (7(2) can be expressed using substitution = x 



kt r°° 

n Jo 



Hz 

g(z) = -— I dxln(l + e- x ) ^ 

" ' Km) +" 2 



Xr 



Furthermore defining ^ = y we obtain, 



XT 



and 



/■OO 

g(y) = -— / tteZna + e-*) 

71" JO 

kt r°° 

9W) = -— Jf d,Mi + e-) r2 + i2 



7/ 2 + X 2 



Therefore, we can write down Eq.44 as 



(3.15) 



(3.16) 



Hv) = g(v) + gtf) 

/ eta; m(l + e 
'0 



7T 



y + y 



y 2 + a; 2 y* 2 + a; 



(3.17) 



This is the central result of this paper which is valid for arbitrary temperature. We now 
consider two different limits (low and high temperature) depending upon the parameter 
value \y\. 
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B. Temperature dependence of Specific heat 



1. Low- Temperature limit: (\y\ 3> 1,KT <C h\z\ 



Partial integration of Eq.46 results in, 



KT f°° tan 1 (-) 

g(y) = / dx- 



7T 



e x + 1 



(3.18) 



When the parameter \y\ ^> 1, corresponding to KT <C h\z\ (which, in turn, implies 
KT <C fkuo), the argument of tan _1 (-) can be expanded as, 



tan 



x 



E 

n=0 



(-1)" fx 



2n + 1 \y 



2n+l 



< 1 



(3.19) 



We then obtain, 

g(y) 



KT 



E 



■i) 



ra+l 



^ ^ ( 2n + l)y 2n+1 Jo e x + l 



x 2n+l 



dx- 



KT 



E 



:-i) 



n+l 



7T ^ (2n + l)y 2n+l 



C(2n + 2)r(2n + 2) 1 



22n+l 



(3.20) 



where ((n) is Reimann-zeta function and T(n) is Gamma function. Now using the above 
results and substituting the values of g(y) and g(y*) in Eq.47, we have 

' z 2n+l _|_ z *(2n+l) ' 



|4n+2 



tt^ (2n+lW n+1 

n=0 v ' 



22n+l / U|4n+2 



*n2?i+1 



5^1 2n+1 ] | z |2fc^2(n-fc)+l + ;z *(2(n-fc)+l)J 

fc=i V k 



(3.21) 



The expression for C v follows immediately, 
d 2 F 



C — —T 

g T 2 



K 



IT 



E 

71=0 



■l) n (2n + 2) (KT 



\An+2 



2n+l 



(z + Z*) 



*\2n+l 



E 

fc=l 



2n + 1 
k 



C(2n + 2)r(2n + 2) 1 - 



22n+l 



2fe f r 2(n-fc)+l _|_ ^*(2(n-fc)+l) 



2 Z 



(3.22) 
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From Eq.52 the specific heat up to three leading order terms are, 




(3.23) 



A close inspection of the above expression clearly reveals that each term in the square bracket 
is a product of two terms. The first one is temperature dependent and of the form (j^^J , 
with n=l,3,5 and originates from traditional thermodynamics. The second one contains 
powers of O^J which has a dynamical origin in k(uj) and is the major content of this work. 
The appearance of dynamical factor is reminiscent of the Kramers theory of activated rate 
process where the rate constant from transition state theory (thermodynamic contribution) 
gets multiplied by this factor in the form f^J where Ub is the frequency of the inverted 
oscillator well at the barrier top. Drawing a hint from this observation, we may therefore 
emphasize that the above expression for C v at low temperature pertains to a small system 
allowing a steady flow of energy from it to the environment in the form of dissipation under 
a quasi-stationary condition. Secondly, we note that the expression for specific heat matches 
with that for the well known form for degenerate fermi system, i.e., C v = AT + BT 3 , (A, 
B being constants), so far as the thermal behaviour of the system is concerned. Our results 
qualitatively agree with the results obtained by Ford and 0'connell|?J and by Hanggi and 
Ingold[8j for a harmonic oscillator in a harmonic bath. The origin of this agreement lies 
on merging of the thermal behaviour of the spin bath and the harmonic bath as T — > as 
discussed earlier. The specific characteristics of the spin system are, however, reflected in 
the details of numerical factors. 



2. High- Temperature limit:(\y\ — > 0,KT 3> Hloq) 

In order to calculate the temperature dependence in the high temperature limit, we 
expand the logarithmic term in Eq.46 as, 

ln(l + e~ x ) = V (-l) n - \e~ x < ll (3.24) 

z — ' n 

n=l 
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Substituting Eq.54 in Eq.46 we obtain, 

TiT'T 1 _ °°_ /*oo 

.'/!//) E< l! " ' ' !h ' 



IT 

KT 

7T 



oo 



y e" 



+ x 2 n 



B-ir 1 



n=l 



Ci(ny)Sin(ny) + -Cos(ny)(n — 2Si(ny)) 



(3.25) 



Making use of Eqs55 in Eq.46 and after formally expanding upto fourth order of y we obtain, 



g(y) = 



KT 



7T 



OO / . x OO 

( 7e - 2)y + y 2 ^(-1)-^ - ( 7 « - 3 ) £(-l) B - V 

n=l ^ ' n=l 



, 2 y_ 

6 



+E(- i r 1 4 + J(i-B-D 



(3.26) 



n=l \ n=l n=l 

where 7 e = 0.577216. After substitution of the values of y, z, z*, the free energy in the high 
temperature regime is given by, 



F(T) 



KT 



7T 



n=l 

3 



/„^ o y I 7e -4/3 



71=1 

OO 

+ E 

n=l 



KT J I 4 

3 



w 



n V ATT / 1 6 

;-l) n_1 fnhwoY J 8 + 7T 

"48" 



n UT 



To 
w 



To 

Wo 

4 



-3 



To 
w 



+ 2-4 



To 
w 



(3.27) 



From Eq.57 the expression for specific heat C v can be calculated as follows; 

d 2 F 



C — —T 

Q T 2 



K 

7Y 



(-1)™ fnfrwoY J (4 — tt) 



n 



E 

71=1 

OO 

+ E 



E 

.71=1 

(~l) n /nW 4 
1 n \KT 

1 



y 1 2 

(Te - 4/3) 



To 
w 



- 2 



7b 
w 



-3 



To 
w 



n=l 



"/'^oV J (8 + 7r) / / "(I V ; j / -;0 



XT 



w 



w 



(3.28) 



Eq.58 shows that the quantities in the square bracket is again a product of an equilibrium 
factor and a dynamical factor that depends on the ratio — . In order to extract out the 
physically relevant form of the above expression we need to consider the thermal saturation 
effect at high temperature since the retention of higher order terms in Eq.58 does not make 
it meaningful above the saturation temperature. This is because the system gets decoupled 
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from the bath above this temperature [T5l [25J [26]. This can be ascertained from the effective 
spectral density and by expressing \tanh (^f) = h — kt , = — ( S z ) where (S z ) is a 



2 £ hu/KT + 1 

measure of the population difference between the two levels of a bath atom. We note that 



as 



(■ 



,tkj/KT 



l) — > 2 the hyperbolic tangent factor tends to vanish as a result of thermal 

C v at high 



saturation of the bath. Keeping therefore only the leading order term of Eq. 
temperature can be expressed as, 

2"! ( ic cs \ / / \ 2 

i i / /v. v \ 

C V = K 




(3.29) 

In order to check this result, we must compare it with the leading order term of specific heat 
of a single spin system in thermal equilibrium. The average energy of a single spin is, 



E = Hu ( S : 



(3.30) 



Substituting (.Sz/ — —\t an h{-^f) in Eq.60 and expanding in inverse powers of T[27J, we 
obtain the leading order term; 



E 



(^o) 2 
ART 



(3.31) 



and the specific heat as, 



C V = K- 



1 / Hou, 



(3.32) 



4 \KT t 

This corresponds to the equilibrium factor in the expression for specific heat at high tem- 
perature in Eq.59 and serves as an important check for our calculation. 



C. Comparison with experiments 

In order to check the validity of our theoretical results we have considered the experi- 
ments dealing with the finite size effects of thin anti-ferromagnetic CoO layers in CoO / S1O2 
multilayer [TU]. Specific heats were evaluated in these experiments as indicators of magnetic 
ordering temperature. It had been demonstrated that specific heat of CoO extracted after 
subtracting the contributions from S1O2 (diamagnetic) from the multilayer with varying 
CoO thickness exhibits a linear variation with temperature upto a maximum followed by a 
decrease. This is shown in Fig. 1. The peak in specific heat curve represents the tempera- 
ture of the maximum change in entropy from a paramagnetic to an anti-ferromagnetic state 
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120 ' 180 240 300 

Temperature 



FIG. 1. (□) and (o) are experimental results of specific heat for 20A and 60A CoO layers, 
respectively [TU]. Solid curves are the results in the low temperature limit and dash curves are 
the results in the high temperature limit of our model for scaled temperature. The ratio 70/wo 
used in our calculation are 1.85 and 2.1 for 20 A and 60 A layers, respectively. 




Diameter (nm) 

200 240 280 320 

Temperature 

FIG. 2. □ (4nm) and o (20nm) are experimental results of specific heat of CoO micrograms [28 . 
■Jo/coq used in our calculations are 1.95 (4nm) and 2.31 (20nm) respectively. [Inset: log-log plot 
of the variation of Neel temperature obtained from the cusp of the theoretical curve for C v as a 
function of grain diameter. Solid line is a two parameter linear fit.] 

for the average CoO grain. Considering the core spins to be coupled to the surface spins of 
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anti-ferromagnetic particles (treated as a spin bath), we have made use of the specific heat 
expressions (14) and (15) for CoO layers with thickness 20A and 60A corresponding to ^ 
ratio 1.85 and 2.1, respectively. The results are plotted in Fig.l by solid and dashed lines, 
respectively. The variation of specific heat with temperature is marked by a divergence or a 
cusp one expects for anti-ferromagnetic transition. For thinner layer the magnetic ordering 
temperature or Neel temperature drops significantly as evident from our analysis presented 
in the figure and observed in experiments. Furthermore the variation of Neel temperature as 
a function of film thickness (d) is depicted in a log- log plot in the inset of Fig. fl. A simple 
two-parameter linear fit corresponds to the relation [2"5] (7& — Tf)/7& — (jM , with T& 
and Tf as the bulk and film Neel temperature is shown. Tf is estimated from the cusp of 
the theoretical curve for C v . The two parameters £o and v are the correlation length at OK 
and critical exponent, respectively. 

The comparison between theory and experiment as discussed above concerns CoO layers. 
The significant reduction in Neel temperature with system size has also been demonstrated 
in a recent experiment [28] on micrograins of CoO. In Fig.f2 we have shown the theoretical 
variation of specific heat with temperature for two different grain sizes, 4nm and 20nm 
and compared with experiments. The specific heat increases monotonically upto a Neel 
temperature which appears as an inflexion of the curve. One observes the decrease of Neel 
temperature with decreasing sample thickness. Finally the variation of Neel temperature 
T/v calculated from theoretical cusp is shown as a function of grain size (diameter) in the 
log-log plot. The critical exponent and correlation length one obtained from the linear fit 
as shown in the inset of Fig. f2. It follows that as one reduces the dimension of the system 
from layer to grain, the correlation length is significantly lowered. In general, the system size 
controls the magnetic properties through the correlation length. The rationale behind the 
choice of higher 70/1^0 ratio for higher thickness of the CoO layers is the following. Since the 
semiconductor band gap is a weak function of layer thickness, ujq remains almost constant 
while it is plausible that electronic conductivity gets hindered by increasing layer thickness. 
This results in higher dissipation and hence larger ratio. Furthermore, a remark on 

accommodating phase transition in finite systems (rather than in the thermodynamic limit) 
in the present analysis is pertinent. The transition temperature is a result of a competition 
between spin alignment and thermal fluctuation. Thermal saturation of the spin bath at 
higher temperatures leads to an effective reduction [91 [HJ US] of system-bath coupling which 
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is again governed by the system size. The fluctuation correction factor therefore plays a 
major role in determining Neel temperature. Although the critical behaviour of correlation 
length is outside the scope of present theory which primarily addresses the problem of 
quantum dissipation in a finite size system, we anticipate this behaviour in the variation 
of Neel temperature (estimated from the theoretical cusp) with system size which can be 
fitted by a power law. The satisfactory agreement between the theoretical and experimental 
results vindicates the validity of our analysis. 

IV. CONCLUSION 

In this work we have examined the behaviour of specific heat of a spin-| system coupled 
to a spin-bath with finite coupling strength. The basis of our analysis is the mapping 
of the spin-spin-bath Hamiltonian into its oscillator counterpart using Holstein-Primakoff 
transformation followed by subsequent diagonalization of the Hamiltonian in its coherent 
state representation. The derived expression for the specific heat is characterized by the 
dynamical correction factors over and above the usual thermal equilibrium factors. By 
tuning the system size it is possible to realize this dynamical contribution to specific heat 
and its consequences. The role of this contribution is elucidated by correlating our theoretical 
analysis to experiments on specific heat of thin metal oxide films. We summarize the main 
conclusions of this study as follows: 

(i) The c-number quantum Langevin equation assumes exactly the traditional classical 
looking form, the essential difference, however, lies on the nature of noise characteristics 
due to the spin-bath. Therefore the susceptibility, the key quantity for calculation of free 
energy associated with the quantum Langevin equation retains also its standard form. The 
underlying universality of the spin bath is therefore reflected in the analysis so far as the 
dynamic behaviour is concerned. 

(ii) The expression for specific heat both in the high and low temperature regions appears 
as a sum of the products of a thermal equilibration factor and a dynamical correction 
factor the later being a function of the ratio of the friction coefficient and the characteristic 
frequency of the system. Up to a leading order the specific heat at low temperature can 
be expressed as C v = A(^ /uo)T + B(^ /ujo)T 3 which reflects the general behaviour of the 
degenerate fermi system. The modification by dissipative contributions due to Ohmic bath 
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is very much similar to that for the bosonic bath. This again highlights the similarities of 
the behaviour of the spin bath and the harmonic bath as T — > 0. 

(iii) The high temperature behaviour of the spin bath is dominated by thermal saturation 
of the two levels of the spin system. Up to a leading order the specific heat below saturation 
temperature can be expressed as C v = C^q/uiq)^. The result corresponds to the qualitative 
behaviour of a single spin-^ system. 
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